Charge doping versus disorder in CeFeAsO: do the in- and out-of-plane dilutions play 

the same role? 
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We provide direct experimental evidence for the identical effect of the in-plane Fei-a;Coa; and 
of the out-of-plane Oi-xFcc chemical dilutions on the itinerant spin-density- wave (SDW) magnetic 
phase in CeFeAsO. Remarkably, the suppression of SDW is not sensitive at all to the different kinds 
of disorder introduced in the two cases. Still, it is clearly shown that the sizeable in-plane disorder 
induced by the Fei_ccCocc substitution is highly effective in suppressing T^. Differently from what is 
observed in CeFeAsOi-^Fx, the ordered magnetic phase of the Ce sublattice is preserved throughout 
the whole phase diagram in CeFei-ccCoxAsO (x < 0.2). An intriguing effect is encountered, whereby 
the magnetic coupling among Ce^^ ions is enhanced by the superconducting phase. 

PACS numbers: 74.70.Xa, 74.25.Dw, 74.25.Ha, 76.75.+i 



The substitution of ~ 10% of 0^~ ions by F~ leads 
to the emergence of high-Tc superconductivity (SC) in 
LaFeAsO [U [2] and in other oxy-pnictides |3J-[£|- Al- 
though such chemical dilution is realized out of the FeAs 
layers, it is known to induce a charge doping in the Fe 
bands. This strategy allows one to reach T^ ^ bh K in 
SmFeAsOi-a^Fa;, namely one of the highest values among 
all the Fe-based SC [THTO]. In spite of the lower T^ val- 
ues, also other kinds of chemical substitutions which still 
give rise to SC in Fe-based materials are widely studied. 
These involve the dilution of transition metal (TM) el- 
ements on the FeAs layers pjj |12J, the most common 
one being probably cobalt p^SHTSl. Whether the chem- 
ical substitution of Fe by other TM elements effectively 
induces a charge doping similarly to the case of Oi-a^F^^ 
has been subject of recent intense efforts on both the 
computational and the experimental sides, particularly 
in the case of Fei-a^Coa^. Contrasting results have been 
reported from density-functional theory (DFT) calcula- 
tions [191 HO] as well as from measurements via x-rays ab- 
sorption [2TJ [22] and photoemission [23 spectroscopies. 
Moreover, the perturbation induced by the Fei_a;TMa^ 
substitution is expected to introduce a high degree of 
in-plane disorder, in turn destabilizing the magnetic itin- 
erant spin-density- wave (SDW) phase [191 El]- This is 
consistent with data for different TM elements in 122, 
where a common scaling for the suppression of SDW is re- 
alized by considering the amount of dilution rather than 
the effective charge doping [11] . Such picture clearly goes 
far beyond the scenario of the weakening of the Fermi- 
surface nesting induced by charge doping [25 . The in- 
terplay of the effects of such disorder on the itinerant 



magnetic phase and SC is even reported to lead to a 
counterintuitive enhancement of T^ with increasing the 
degree of disorder in the underdoped region of the phase 
diagram [26 . 

In this Letter we provide direct evidence for the equiva- 
lence of the in-plane Fei-a^Coa^ and out-of-plane Oi-a^F^; 
chemical dilutions in CeFeAsO as far as the magnetic 
properties are concerned. Measurements of muon-spin 
spectroscopy (/i+SR) show how both these substitutions 
lead to a quantitatively identical suppression of the SDW 
phase as a function of x. Two different regimes are 
detected and a short-range magnetic phase is found to 
nanoscopically coexist with SC within a narrow range of 
X values. Remarkably, the SDW transition temperature 
Tn is not affected at all by the degree of the in-plane 
disorder. However, the in-plane disorder itself strongly 
affects the properties of SC, T^ being sizeably suppressed 
by the Fei-a^Co^^ dilution with respect to what is real- 
ized via the Oi-^^x substitution. The magnetic order- 
ing temperature T§^ for the Ce sublattice is preserved 
throughout the whole phase diagram in CeFei-a^Coa^AsO 
independently from the SDW phase, at variance with 
what is reported for CeFeAsO i-a^F^^ [6l|10]. Finally, an 
intriguing correlation among the dome-like features of 
the phase diagram for both T^ and T^^ is described sug- 
gesting an enhancement of the magnetic coupling among 
Ce^+ ions mediated by the SC phase. 

Powder samples of CeFei-a^Co^^AsO (0 < x < 0.2) 
were synthesized as described in detail in the supplemen- 
tal material (SM). The characterization of the samples 
was performed using electrical transport [27 and dc mag- 
netometry, while a detailed investigation of Tn was per- 
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FIG. 1: (Color online) ZF /x^-depolarization curves for four 
CeFei-ccCoccAsO samples (0 < x < 0.07) at representative T 
values. Continuous lines are best fits to experimental data 
according to the analysis described in the SM. Data for x = 
are taken from Ref. 6 . 



FIG. 2: (Color online) Electronic phase diagram for 
CeFei-ccCoccAsO. Mpe and Mce refer to magnetic phases of 
Fe and Ce, respectively. The dashed lines are empirical guides 
for the eyes. Dashed lines on the magnetic side of the phase 
diagram are discussed in more detail in Fig. [S] 



formed by means of /i+SR [28]. The results of the char- 
acterizations as well as the precise information about the 
definitions of the critical temperatures Tn, T^ and T^^ 
are reported in the SM. 

Sets of experimental /i+ spin-depolarization curves 
AT{t) in zero-magnetic field (ZF) at different tempera- 
tures (T) are shown in Fig.[l]for several CeFe i-a^Co^^ AsO 
samples {0 < x < 0.07). Data are analyzed as described 
in the SM. Two inequivalent crystallographic sites for /i+ 
are observed in agreement with previous results on oxy- 
pnictides [29H3T] . Because of the dipolar interaction be- 



tween /i+ and the antiferromagnetic FeAs layers, /i+ are 
very sensitive to magnetic ordering even when the mag- 
netic coherence length ^^ is of the order of few nanome- 
ters. Typically, the precessions of /i+ around the local 
magnetic fields are rapidly overdamped when ^m ^ 10 
lattice constants [28l [S^. This is exactly the case of 
CeFei_a^Coa;AsO for x > 0.035, as shown in Fig. [l] a re- 
sult indicating that the long-range magnetic order (LRO) 
realized at x = 0, where ^m ^ 1 nm, is gradually driven 
to a short-range magnetic order (SRO), where ^m "^ 1 
nm for higher x concentrations. 

The full electronic phase diagram of CeFe i-a^Co^^ AsO 
is presented in Fig. [2] By comparing these results with 
what is reported for CeFeAsOi-a^Fa^ |^, it is immediately 
clear that the SC properties are affected differently by the 
two chemical dilutions. In particular, the typical values 
Tc :^ 10 K for CeFei-a^Coa^AsO at optimal doping are 
much lower than what is observed for CeFeAsOi-a^F^^, 
namely Tc ::^ 35 K [1 [lOl US]. This suggests that the 
Fei_a;Coaj dilution is generally detrimental to SC since it 
introduces a sizeable degree of disorder directly on the 
active FeAs layers. However, the region of coexistence 
between SDW and SC phases is similar to what was re- 
ported for SmFeAsOi-a^Fa; and CeFeAsOi-a^F^^ ^iS^. In 



particular, an extremely narrow coexistence region for 
the two electronic ground states is observed for both the 
in-plane and out-of-plane substitutions in the 1111 fam- 
ily. Under these conditions, 100% of the sample volume 
is magnetic while, at the same time, a sizeable bulk re- 
gion displays SC (see SM). From these arguments, a seg- 
regation of the two phases at the nanoscopic level was 
deduced in SmFeAsOi-a^F^^ and CeFeAsOi-a^Fa^ and the 
same scenario holds also in the case of Fei-a^Coa^ dilu- 
tion. This scenario is completely different from what is 
detected in La-based 1111 materials [2][33l|34]. It should 
be stressed that such nanoscopic coexistence is detected 
in CeFeAsOi-a^Fa; and CeFei-a^Coa^AsO only when the 
magnetic phase is short-ranged. Remarkably, this is qual- 
itatively different from what is obtained in the electron- 
doped 122 materials like Ba(Fei_a;Coa^)2As2, where the 
x-range of coexistence is typically much wider [TTJ |35] , 
superconductivity being reported to emerge still in the 
presence of clear coherent oscillations for the SDW phase 
[36]. The appearance of superconductivity even in the 
presence of a LRO magnetic phase is possibly a finger- 
print of the 122 family, at variance with what is gener- 
ally observed in 1111 materials. This can be observed by 
comparing the results presented in Ref. [36 with what 
is reported in the case of hole-doping realized by means 
of the Bai_a^Ka^ substitution [37]. 

With the aim to precisely compare the effects of the 
in-plane and out-of-plane dilutions in CeFe AsO, the val- 
ues of Tn versus x in CeFe i-a^Co^^ AsO are reported 
in the main panel of Fig. |3) together with data for 
CeFe AsO i-a^Fa^, the latter being taken from Ref. [6j. 
It should be stressed that x concentrations for the latter 
materials should be considered as real after an a poste- 
riori determination by means of nuclear magnetic reso- 
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FIG. 3: (Color online) Main panel: Tn versus x for Co- 
and F-doped CeFeAsO (the latter data are taken from Ref. 
[6]). Close and open symbols refer to estimates from /x^SR 
and resistivity, respectively. The shaded region indicates the 
crossover between the LRO and SRO magnetic regimes. The 
dashed lines are empirical guides for the eyes. Inset: compar- 
ison between ZF /i+ -depolarization curves at representative 
T for samples in the SC region of the phase diagram for both 
CeFei-xCoccAsO (x = 0.12), panel (a), and CeFeAsO i-^Fa; 
{x = 0.08), panel (b). Continuous lines are best fits to exper- 
imental data according to the analysis described in SM. 



nance [6]. Remarkably, results for both the families of 
compounds display a common behavior. Magnetism is 
indeed driven only by the amount of dilution regardless 
of the possible different modifications both of the de- 
tails of the band structure or of the introduced degree of 
in-plane disorder. Two qualitatively different trends for 
the suppression of Tn with increasing x can be detected, 
which are associated with LRO and SRO (denoted as 
Ilro and IIsro in Fig. [3J respectively). In particular, a 
linear dependence of Tn on x common to both F- and Co- 
diluted compounds is encountered in the intermediate- 
doping region close to the appearance of SC. The results 
are similar to what is realized in Fei+^Se^^Tei-a^, where 
the region II is associated with the emergence of a glassy- 
like magnetism [38l |39] . Remarkably, the phase diagram 
is also highly reminiscent of hole-doped cuprates, such 
as La2-xSrajCu04 and Yi_a^Caa;Ba2Cu306+ty5 where the 
glassy magnetism is also generally reported to be fully 
suppressed for x ~ 0.1 [40l|41]. However, the coherent 
oscillations of the /i+SR signal are clearly detected also 
well-inside region II both for Fei+^ySea^Tei-a^ and for the 
cuprates due to clusters whose magnetic moment freezes 
at low T [39, ,41j . This inhomogeneous magnetic order is 
often referred to as cluster spin glass. 

In spite of the above similarities, clear qualitative 
differences among the effects of in-plane and out-of- 
plane disorders are enlightened concerning the features 
of the magnetic phase of the Ce sublattice. As reported 
in Refs. [6^ and [lO^, T^^ is strongly suppressed in 



CeFeAsO i_ a; Fa; as x increases and eventually vanishes 
in correspondence with the full suppression of the SDW 
phase. A similar phenomenology was recently reported 
in CeFci-a^Rua^AsO even though a much weaker depen- 
dence on X was detected in this case \42^. This would 
suggest that the two magnetic phases are intimately con- 
nected and that the ordering of the Ce sublattice is in- 
duced by the molecular field generated by the SDW [lOj. 
The current results clarify that this is not really the case, 
since data for CeFei-a^Co^^AsO unambiguously show that 
the magnetic phase of the Ce sublattice is still present at 
X = 0.2, although the SDW is fully suppressed already 
at X = 0.1 (see Figs. |2]and[3|. Similar results were re- 
cently reported also in the case of isovalent dilution in 
CeFeAsi-a^Pa^O [43]. This observation for samples well 
inside the SC phase is in agreement with what is reported 
in Ref. [16 and it is here reinforced by /i+SR data (see 
Fig. [3J inset). One clearly detects coherent oscillations 
in the case of CeFeo.88Coo.i2AsO indicative of a long- 
range ordered phase, even though the transverse damp- 
ing is quite severe {Xj^ c^ 4 /is), see Fig. pfa). On the 
other hand, in CeFeAsO0.92F0.08 (nominal composition) 
the ordering of Ce^+ ions can be detected from a qualita- 
tive change in the damping of the transverse relaxation, 
which from Lorentzian turns into a markedly Gaussian 
shape at the lowest T value, see Fig. IsFb). No clear signs 
of coherent oscillations can be distinguished in this case, 
again confirming that the Oi-^^F^; dilution strongly af- 
fects the magnetic phase of the Ce^+ ions. 

These results can be interpreted in the light of two 
different mechanisms contributing to the magnetic cou- 
pling among Ce+ ions, one of direct nature (e.g., superex- 
change via 0^~) and the other of indirect RKKY-like 
nature, via carriers belonging to the FeAs layers |44| l45]. 
The strong sensitivity of T^^ to the local perturbation 
realized by the Oi-a^Fa^ dilution shows that the former 
mechanism has the strongest impact to the overall mag- 
netic coupling. Naively, this is favored by the spatial 
extension of the Ce^+ electronic orbit als, typically larger 
than those of other rare-earth ions. The sizeable sen- 
sitivity of the direct term to the increasing F content 
completely hinders the modification of the RKKY inter- 
action realized by the variation of the charge density on 
the FeAs bands and of the relative band structure. On 
the other hand, the in-plane Fei-a^Co^; dilution does not 
affect the CeO layers locally and the main modification 
is expected to involve the indirect term. Moreover, the 
change of the extent of the f — d hybridization among 
Ce^+ ions and FeAs bands, triggered by the reduction 
of the volume of the crystallographic cell with increas- 
ing X [16], should also be considered. In this respect, it 
should be noticed that the gradual increase of T§^ upon 
increasing x in CeFei-a^Coa^AsO (see the dashed straight 
line in Fig. [4]) is in qualitative agreement with the simi- 
lar trend observed after /i+SR measurements in undoped 
CeFeAsO under applied hydrostatic pressure [30 . These 
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FIG. 4: (Color online) Main panel: Tc values as reported after 
measurements of dc magnetometry. T^ ^ is also plotted (on a 
different scale) for two different values of magnetic field. The 
continuous and dashed lines are empirical guides for the eyes. 



results deserve a more detailed computational investiga- 
tion in order to precisely unveil the details of such com- 
plex interplay of subtle many-body effects. 

Quite surprisingly, the magnetic coupling among Ce^+ 
ions is anomalously enhanced in the presence of SC. 
This is clearly evidenced in Fig. [4j where a dome-like 
behaviour for T^^ versus x in correspondence with the 
emergence of the SC dome is superimposed to the overall 
linear increase described above. T^^ deduced by means of 
dc magnetometry is known to be strongly suppressed as 
the value of the measuring field increases for typical val- 
ues H ^ 10 kOe (see SM). At first sight, the enhancement 
of Tj^® could be associated with the shielding of the ex- 
ternal field by the SC-induced supercurrents. This effect 
would apparently increase the critical temperature itself. 
However, this scenario is ruled out by the tiny discrep- 
ancy observed among measurements at H = 20 Oe and 
10 kOe, thus implying that the applied H field value can- 
not be at the origin of the observed T^^ suppression. The 
explanation possibly relies on the modification of the in- 
direct RKKY coupling through the involvement of FeAs 
bands into SC. Still, the results are even more intriguing 
in view of the possible occurrence of subtle many-body 
effects like the Kondo screening and its interplay with 
SC. This issue has been discussed in the literature from 
a theoretical point of view, showing how a competition 
between these two mechanisms is at work in oxy-pnictides 
[46^. Our results display a common enhancement of Ce 
magnetism and SC and are not in contrast with such 
theoretical scenario. 

Summarizing, in this Letter we discussed the effects of 
the in-plane and out-of-plane chemical substitutions in 
CeFeAsO realized by Fei-^^Coa^ and Oi-a^Fa^ dilutions, 
respectively. Remarkably, /i+SR enlightened that the 
spin-density-wave phase is suppressed in a quantitatively 



identical fashion in both cases, showing that the degree 
of the in-plane disorder is not playing any substantial 
role. At variance with the behavior for the spin-density- 
wave, such in-plane disorder strongly affects the super- 
conducting properties resulting in a much more effective 
suppression of Tc for the Fei_a;Coa; rather than for the 
Oi-a^Fa^ dilution. The unusual behavior of Ce magnetism 
was discussed in the light of the different kinds of disor- 
der introduced by the two considered chemical substitu- 
tions. Our surprising results clearly deserve a more de- 
tailed computational investigation concerning the precise 
role of indirect RKKY-like coupling among the localized 
Ce+ magnetic moments, also in relation with the anoma- 
lous enhancement of the ordering temperature T^^ in the 
presence of superconductivity. 
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SUPPLEMENTAL MATERIAL 

Synthesis of samples 

Polycrystalline samples of CeFei-a^Coa^AsO {0 < x < 
0.2) have been prepared by a two-step solid-state reac- 
tion similarly to what is described in Ref. [47 . In the 
first step, CeAs was prepared from Ce slugs (Chempur, 
99.9%) and As lumps (Alfa Aesar, 99.999%) reacting a 
stoichiometric ratio in an evacuated quartz tube placed 
in a two-zone furnace. In the second step, we used 
the resulting CeAs and mixed it with Fe (Alfa Aesar, 
99.998%), FesOa (Chempur, 99.999%), and Co (Good- 
fellow, 99.99%) in a stoichiometric ratio. All starting 
materials were homogenized by grinding in a ball mill. 
The resulting powders were pressed into pellets under Ar 
atmosphere using a pressure of 20 kN, and subsequently 
annealed in an evacuated quartz tube in a two-step an- 
neahng process at 940° C for 8 h and at 1150°C for 48 
h. 

In order to confirm the single-phase character of the 
poly crystals, powder x-ray diffraction was performed on 
a Huber Guinier camera (Co K^ radiation). The samples 
were either pure or contained small amounts of CeAs, 




FIG. 5: (Color online) p versus T for CeFei-ccCoxAsO sam- 
ples in ZF conditions. Data for non-SC and SC samples are 
plotted in the upper and lower panels, respectively. Upper in- 
set: enlargement of data for non-SC samples displaying kinks 
at Tn^ (see arrow - different scales are used for the two sam- 
ples). Lower inset: enlargement of data for SC samples in the 
low-T region. 



C03O4, Fe304 and/or FeAs. The microstructure and 
the composition were examined by scanning electron mi- 
croscopy (XL30 Philipps, IN400) equipped with an elec- 
tron microprobe analyzer for semi-quantitative elemental 
analysis using the wavelength dispersive x-ray mode. The 
analysis showed a good agreement between the nominal 
and the measured Co contents. Accordingly, the nominal 
content is used to label the samples. 



Characterization of samples 

Electrical transport 

Measurements of electrical resistivity (p) were per- 
formed in conditions of zero-magnetic field (ZF) as a 
function of temperature (T) by means of a standard four- 
probe setup. The T-dependence of p for CeFei-a^Coa^AsO 
samples (0 < x < 0.15) is shown in Fig. [s] The curves 
are normalized with respect to the relative prt values at 
room T (RT) for a better visualization of data. In the 
case of CeFeAsO, the pronounced maximum at around 
150 K and the inflection point at slightly lower T can 
be related to the critical temperatures for the structural 
and spin-density- wave (SDW) transitions, Tg and Tn re- 
spectively [TOl [27] . These anomalies can be discerned up 
to X = 0.03 even if a sizeable and progressive broadening 
prevents one from a precise definition of Tg and Tn for 
X > 0.01. However, a progressive suppression of both 
Tg and Tn with increasing x can be inferred as a clear 
trend, in agreement with data reported in the literature 
[16] . The clear kinks in p versus T curves at T c^ 4 K are 
associated to the magnetic ordering temperature T^^ for 




FIG. 6: (Color online) dc susceptibility for CeFei-ccCo^AsO 
SC samples. Measurements were performed at if = 20 Oe 
in ZF-cooled conditions. The clear kinks at T 2^ 4 K are 
associated with T§^ (see also Fig.m). 



the Ce sublattice (see the upper inset of Fig. [5|. A fur- 
ther increase in the level of Co-doping (0.07 < x < 0.12) 
leads to the appearance of superconductivity (SC) with 
the highest value for the critical temperature Tc being re- 
alized for the optimal-doping value x !:^ 0.1 (see the inset 
in the lower panel of Fig. |5| while SC is fully suppressed 
by increasing the Co doping to values x > 0.15. Our data 
are in good agreement with early reports in the literature 



[T4HT6]. 



dc magnetometry 



The magnetic characterization of the samples was per- 
formed by means of a MPMS-XL5 SQUID dc magne- 
tometer (Quantum Design). Results for the SC materi- 
als are presented in Fig. |6] where a dome-like behaviour 
for the critical temperature T^ for the SC phase is evi- 
dent. The maximum shielding fraction can be roughly 
estimated as 30%, in good agreement with previous re- 
ports on F-doped CeFeAsO [6l |32] even if systematically 
lower values are found in the current case (no demagne- 
tization effects from the geometrical shape of the grains 
were considered at all). However, the picture is very use- 
ful in order to check that the SC phase is not filamentary 
in the considered materials even if a precise estimate of 
the SC fraction can not be obtained from the current 
data. 

The value of T§^ can be clearly detected by means of 
dc magnetometry since magnetization M displays sharp 
peaks in correspondence to those critical temperatures 
(see figs. [6]and[7|). A strong dependence of such quan- 
tity on the external magnetic field is detected as shown 
in the main panel of Fig. [Tj However, it should be re- 
marked that in the limit JY < 10 kOe the results display 
an almost negligible dependence of T§^ on H (see the in- 
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FIG. 7: (Color online) Main panel: dc susceptibility for Ce- 
FeAsO in the high- if regime. The arrows indicate T§^. Inset: 
dc susceptibility for CeFeo.8Coo.2AsO in the low- if regime. 



set of Fig.!?]). The error bars associated to results for SC 
samples at low- if values are by far too big if compared 
to what is obtained for non-SC materials (see Fig. [6|. 
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FIG. 8: (Color online) Upper panel: T-dependence of the 
internal field at the /x^ site 5^^ for x = and x = 0.03 
(both F- and Co-dilutions are considered). The continuous 
lines are best-fits to data according to eq. ([2]). Lower panel: 
T-dependence of the magnetic volume fraction Vm (T) for the 
CeFei-ccCoccAsO samples (x < 0.07). The dashed lines are 
best-fits to data according to eq. (l3|. 



Muon spin spectroscopy 

In a /i+SR experiment, a spin-polarized beam of pos- 
itive muons /i+ is implanted into the investigated sam- 
ple. Coulomb scattering processes typically lead to the 
thermalization of /i+ at interstitial crystallographic sites 
without the loss of spin polarization. After a mean life 
time r^+ c^ 2.17/is, every /i+ decays emitting a positron 
e+ with linear momentum preferentially parallel to the 
direction of the spin of the /i+ at the instant of the decay. 
As a consequence, a spatially-resolved detection of e+ as 
a function of time (t) allows one to probe magnetism 
on a local scale [28]. The typical output of a ZF-/i+SR 
experiment for magnetic materials is the depolarization 
function 
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describing the t-dependence of the spin polarization of 
/i+ at T. Here, Aq is an instrument-dependent parameter 
corresponding to the condition of full spin polarization. 
In the paramagnetic phase of the investigated material, 
/i+ are subject to nuclear magnetism leading to a slow 
Gaussian damping quantified by (Jn- In the presence of 
a magnetic phase, a fraction Vm{T) of /i+ probe a static 
local magnetic field and this quantity reflects the mag- 
netic volume fraction of the sample at T, accordingly. 
The labels Tr and L in eq. (fTl) refer to /i+ probing local 



magnetic fields B^ in transverse or longitudinal direc- 
tions with respect to the initial spin polarization, respec- 
tively. F{t) oscillating functions represent the precession 
of /i+ around the local field 5^ while D^^'^{t) quantify 
the damping of the signal. D^^{t) is typically associated 
with the static distribution of local magnetic fields while 
D^(t) probes 1/Ti-like dynamical processes. The index 
i in eq. (fTl) allows for the possible presence of inequiv- 
alent crystallographic sites for /i+ whose population is 
controlled by the parameters ("i such that ^^ = 1- 

Measurements of /i+SR were performed at S/iS (Paul 
Scherrer Institute, Switzerland) on both the GPS and 
Dolly low-background spectrometers at different temper- 
atures (T) and in ZF conditions. The presence of two in- 
equivalent crystallographic sites for /i+ is envisaged (i = 
1,2) in agreement with previous results on oxy-pnictides 
[29H3T] . In the current case of CeFei_a^Goa;AsO, one finds 
C1/C2 — 4 for the occupation probabilities of the sites Q 
independently from the actual x value. The fitting results 
show that the magnetic volume fraction Vm{T) associated 
with the SDW phase is maintained at around 100% in the 
low-T saturation limit independently from x (see Fig. [sj 
lower panel). /i+ implanted at site 1 probe long-range 
magnetism in CeFeAsO as clearly seen by the long-lived 
coherent oscillation in the Tr component. In particular, 
the function Fiit) = cos (j^Bfj^^t) describes well the data 
while Fi{t) = Jo (7^ 5^1 1) is needed in order to repro- 
duce the experimental trend for x = 0.03 (Jo standing 
for a zeroth-order first-kind Bessel function and 7^ being 
the gyromagnetic ratio for /i+). This is a well-known re- 
sult for oxy-pnictide materials and it hints at the gradual 
modification of the SDW phase from commensurability 



to incommensurability with the underlying lattice upon 
increasing the charge doping while at the same time re- 
ducing the nesting of the Fermi surface 23 HH]- The 
results for B^.^ as a function of T are displayed in the 
upper panel of Fig. [S] In the same figure, it is shown 
how the phenomenological function 



B^,{T) = Blf 



1- ^ 



T 



C'cw \ 

T-e) 



(2) 



first proposed in Ref. ^ correctly reproduce the exper- 
imental data also in the case of the Co-doped x = 0.03 
sample. In particular, the first term on the right side 
of eq. Q accounts for a double-exponent power-law-like 
behaviour of the internal field due to the magnetic order- 
ing of Fe while the second term accounts for a feedback 
paramagnetic polarization of the Ce sublattice induced 
by the molecular field of Fe and governed by the Curie- 
Weiss constant Ccw and the Curie- Weiss temperature 9 
[29 . In the fitting procedure for all the three samples, 
the values a = 2.1 and /3 = 0.2 have been kept fixed 
in order to reduce the number of free parameters (the 
values have been extracted from what is reported in the 
literature [29]). 

The gradual disordering of the magnetic phase leads to 
the disappearance of oscillations in the Tr component for 
X > 0.05, Fi{t) being overdamped by a sizeable distri- 
bution of static local magnetic fields. This corresponds 
to Fi{t) = 1 for both x = 0.05 and 0.07 samples with 
transversal dampings described as Dj^{t) = exp {—Xj^t). 
At the lowest investigated temperature (T = 5 K), one 
has Xj^ c^ 45 jas and 25 /is for x = 0.05 and 0.07, 
respectively. Such values are quite typical in the re- 
gion of coexistence between magnetism and SC for oxy- 
pnictides ^i32j. On the other hand, /i+ implanted at 
site 2 probe sizeable distributions of local magnetic fields 
for all the samples and F2{t) = 1, accordingly, with 



D^'{t) = exp(-Aj^t). At T 



5 K, one has Aj^ c^ 25 



/iS, 20 /is, 15 /is and 8 /as for x = 0, 0.03, 0.05 and 0.07, 
respectively. 

Tn is precisely quantified from the actual T- 
dependence of the magnetic volume fraction Vjji{T) after 
a fitting procedure according to the error-function-like 
expression [611301132] 



Kn(T) = i erfc 






(3) 



where the complementary error function erfc(x) is defined 



as 



erfc(x) 



2 r+^ 



'dt. 



(4) 



In particular, Tn turns out to be defined as the T value 
corresponding to 50% of the magnetic volume fraction. 



* E-mail: g.prando@ifw-dresden.del 
[1] Y. Kamihara, T. Watanabe, M. Hirano, H. Hosono, 

Journ. Am. Chem. Soc. 130, 3296 (2008) 
[2] H. Luetkens, H.-H. Klauss, M. Kraken, F. J. Litterst, T. 

Dellmann, R. Klingeler, C. Hess, R. Khasanov, A. Am- 

ato, C. Baines, M. Kosmala, O. J. Schumann, M. Braden, 

J. Hamann-Borrero, N. Leps, A.Kondrat, G. Behr, J. 

Werner, B. Biichner, Nature Mat. 8, 305 (2009) 
[3] J. Zhao, Q. Huang, C. de la Cruz, S. Li, J. W. Lynn, Y. 

Chen, M. A. Green, G. F. Chen, G. Li, Z. Li, J. L. Luo, 

N. L. Wang, R Dai, Nature Mat. 7, 953 (2008) 
[4] A. J. Drew, C. Niedermayer, P. J. Baker, F. L. Pratt, S. 

J. Blundell, T. Lancaster, R. H. Liu, G. Wu, X. H. Chen, 

I. Watanabe, V. K. Malik, A. Dubroka, M. Rossle, K. W. 

Kim, C. Baines, C. Bernhard, Nature Mat. 8, 310 (2009) 
[5] S. Sanna, R. De Renzi, G. Lamura, C. Ferdeghini, A. 

Palenzona, M. Putti, M. Tropeano, and T. Shiroka, Phys. 

Rev. B 80, 052503 (2009) 
[6] T. Shiroka, G. Lamura, S. Sanna, G. Prando, R. De 

Renzi, M. Tropeano, M. R. Cimberle, A. Martinelli, C. 

Bernini, A. Palenzona, R. Fittipaldi, A. Vecchione, P. 

Carretta, A. S. Siri, C. Ferdeghini, and M. Putti, Phys. 

Rev. B 84, 195123 (2011) 
[7] Z.-A. Ren, W. Lu, J. Yang, W. Yi, X.-L. Shen, Z.-C. Li, 

G.-C. Che, X.-L. Dong, L.-L. Sun, F. Zhou, and Z.-X. 

Zhao, Chin. Phys. Lett. 25, 2215 (2008) 
[8] D. C. Johnston, Adv. Phys. 59, 803 (2010) 
[9] G. Prando, A. Lascialfari, A. Rigamonti, L. Romano, S. 

Sanna, M. Putti, and M. Tropeano, Phys. Rev. B 84, 

064507 (2011) 
[10] H. Maeter, J. E. Hamann-Borrero, T. Goltz, J. Spehling, 

A. Kwadrin, A. Kondrat, L. Veyrat, G. Lang, H.-J. Grafe, 

C. Hess, G. Behr, B. Biichner, H. Luetkens, C. Baines, 
A. Amato, N. Leps, R. Klingeler, R. Feyerherm, D. Ar- 
gyriou, and H.-H. Klauss, arXiv:1210.6959 (2012) 

[11] N. Ni, A. Thaler, J. Q. Yan, A. Kracher, E. Colombier, 
S. L. Bud'ko, P. C. Canfield, and S. T. Hannahs, Phys. 
Rev. B 82, 024519 (2010) 

[12] M. A. Tanatar, N. Ni, A. Thaler, S. L. Bud'ko, P. C. Can- 
field, and R. Prozorov, Phys. Rev. B 84, 014519 (2011) 

[13] A. S. Sefat, A. Huq, M. A. McGuire, R. Jin, B. C. Sales, 

D. Mandrus, L. M. D. Cranswick, P. W. Stephens, and 
K. H. Stone, Phys. Rev. B 78, 104505 (2008) 

[14] J. Prakash, S. J. Singh, S. Patnaik, and A. K. Ganguli, 

Solid State Commun. 149, 181 (2009) 
[15] L.-D. Zhao, D. Berardan, C. Byl, L. Pinsard-Gaudart, 

and N. Dragoe, J. Phys.: Condens. Matter 22, 115701 

(2010) 
[16] T. Shang, L. Yang, Y. Chen, N. Cornell, F. Ronning, J. 

L. Zhang, L. Jiao, Y. H. Chen, J. Chen, A. Howard, J. 

Dai, J. D. Thompson, A. Zakhidov, M. B. Salamon, and 

H. Q. Yuan, Phys. Rev. B 87, 075148 (2013) 
[17] A. S. Sefat, R. Jin, M. A. McGuire, B. C. Sales, D. J. 

Singh, and D. Mandrus, Phys. Rev. Lett. 101, 117004 

(2008) 
[18] A. Yamamoto, J. Jaroszynski, C. Tarantini, L. Balicas, 

J. Jiang, A. Gurevich, D. C. Larbalestier, R. Jin, A. S. 

Sefat, M. A. McGuire, B. C. Sales, D. K. Christen, D. 

Mandrus, Appl. Phys. Lett. 94, 062511 (2009) 
[19] H. Wadati, L Elfimov, and G. A. Sawatzky, Phys. Rev. 

Lett. 105, 157004 (2010) 



[20] T. Berlijn, C.-H. Lin, W. Garber, and W. Ku, Phys. Rev. 
Lett 108, 207003 (2012) 

[21] E. M. Bittar, C. Adriano, T. M. Garitezi, P. F. S. Rosa, L. 
Mendonca-Ferreira, F. Garcia, G. de M. Azevedo, P. G. 
Pagliuso, and E. Granado, Phys. Rev. Lett. 107, 267402 
(2011) 

[22] M. Merz, F. Filers, T. Wolf, P. Nagel, H. v. Lohneysen, 
and S. Schuppler, Phys. Rev. B 86, 104503 (2012) 

[23] G. Levy, R. Sutarto, D. Chevrier, T. Regier, R. Blyth, 
J. Geek, S. Wurmehl, L. Harnagea, H. Wadati, T. Mi- 
zokawa, L S. Elfimov, A. Damascelli, and G. A. Sawatzky, 
Phys. Rev. Lett. 109, 077001 (2012) 

[24] P. Bonfa, P. Carretta, S. Sanna, G. Lamura, G. Prando, 
A. Martinelli, A. Palenzona, M. Tropeano, M. Putti, and 
R. De Renzi, Phys. Rev. B 85, 054518 (2012) 

[25] C. Liu, G. D. Samolyuk, Y. Lee, N. Ni, T. Kondo, A. 
F. Santander-Syro, S. L. Bud'ko, J. L. McChesney, E. 
Rotenberg, T. Valla, A. V. Fedorov, P. C. Canfield, B. N. 
Harmon, and A. Kaminski, Phys. Rev. Lett. 101, 177005 
(2008) 

[26] R. M. Fernandes, M. G. Vavilov, and A. V. Chubukov, 
Phys. Rev. B 85, 140512 (2012) 

[27] C. Hess, A. Kondrat, A. Narduzzo, J. E. Hamann- 
Borrero, R. Klingeler, J. Werner, G. Behr, and B. 
Biichner, Europhys. Lett. 87, 17005 (2009) 

[28] A. Yaouanc and P. Dalmas de Reotier, Muon Spin Ro- 
tation, Relaxation, and Resonance: Applications to Con- 
densed Matter, Oxford University Press, Oxford (2011) 

[29] H. Maeter, H. Luetkens, Y. G. Pashkevich, A. Kwadrin, 
R. Khasanov, A. Amato, A. A. Gusev, K. V. Lamonova, 
D. A. Chervinskii, R. Klingeler, C. Hess, G. Behr, B. 
Biichner, and H.-H. Klauss, Phys. Rev. B 80, 094524 
(2009) 

[30] R. De Renzi, P. Bonfa, M. Mazzani, S. Sanna, G. Prando, 
P. Carretta, R. Khasanov, A. Amato, H. Luetkens, M. 
Bendele, F. Bernardini, S. Massidda, A. Palenzona, M. 
Tropeano, and M. Vignolo, Supercond. Sci. Technol. 25, 
084009 (2012) 

[31] G. Prando, P. Bonfa, G. Profeta, R. Khasanov, F. 
Bernardini, M. Mazzani, E. M. Briining, A. Pal, V. P. 
S. Awana, H.-J. Grafe, B. Biichner, R. De Renzi, P. Car- 
retta, and S. Sanna, Phys. Rev. B 87, 064401 (2013) 

[32] S. Sanna, R. De Renzi, T. Shiroka, G. Lamura, G. 
Prando, P. Carretta, M. Putti, A. Martinelh, M. R. Cim- 
berle, M. Tropeano, and A. Palenzona, Phys. Rev. B 82, 
060508 (2010) 

[33] R. Khasanov, S. Sanna, G. Prando, Z. Shermadini, M. 
Bendele, A. Amato, P. Carretta, R. De Renzi, J. Karpin- 
ski, S. Katrych, H. Luetkens, and N. D. Zhigadlo, Phys. 
Rev. B 84, 100501 (2011) 

[34] G. Prando, S. Sanna, G. Lamura, T. Shiroka, M. Tro- 



peano, A. Palenzona, H.-J. Grafe, B. Biichner, P. Car- 
retta, and R. De Renzi, Phys. Stat. Sol. B 250, 599 
(2013) 

[35] R. M. Fernandes, D. K. Pratt, W. Tian, J. Zarestky, A. 
Kreyssig, S. Nandi, M. G. Kim, A. Thaler, N. Ni, P. 
C. Canfield, R. J. McQueeney, J. Schmalian, and A. L 
Goldman, Phys. Rev. B 81, 140501 (2010) 

[36] C. Bernhard, C. N. Wang, L. Nuccio, L. Schulz, O. Za- 
harko, J. Larsen, C. Aristizabal, M. Willis, A. J. Drew, 
G. D. Varma, T. Wolf, and C. Niedermayer, Phys. Rev. 
B 86, 184509 (2012) 

[37] E. Wiesenmayer, H. Luetkens, G. Pascua, R. Khasanov, 
A. Amato, H. Potts, B. Banusch, H.-H. Klauss, and D. 
Johrendt, Phys. Rev. Lett. 107, 237001 (2011) 

[38] N. Katayama, S. Ji, D. Louca, S.-H. Lee, M. Fujita, T. J. 
Sato, J. S. Wen, Z. J. Xu, G. D. Gu, G. Xu, Z. W. Lin, 
M. Enoki, S. Chang, K. Yamada, and J. M. Tranquada, 
Journ. Phys. Soc. Jpn. 79, 113702 (2010) 

[39] G. Lamura, T. Shiroka, P. Bonfa, S. Sanna, F. Bernar- 
dini, R. De Renzi, R. Viennois, E. Giannini, A. Piriou, N. 
Emery, M. Cimberle, and M. Putti, J. Phys.: Condens. 
Matter 25, 156004 (2013) 

[40] M.-H. Julien, Physica B 329-333, 693 (2003) and refer- 
ences therein 

[41] S. Sanna, F. Coneri, A. Rigoldi, G. Concas, S. R. Giblin, 
and R. De Renzi, Phys. Rev. B 82, 100503 (2010) 

[42] C. Wang, H. Jiang, Y. Luo, C. Feng, W. Li, Z. Xu, G. 
Cao, J. H. Kim, and S. Dou, Europhys. Lett. 99, 57009 
(2012) 

[43] A. Jesche, T. Forster, J. Spehling, M. Nicklas, M. de 
Souza, R. Gumeniuk, H. Luetkens, T. Goltz, C. Krellner, 
M. Lang, J. Sichelschmidt, H.-H. Klauss, and C. Geibel, 
Phys. Rev. B 86, 020501 (2012) 

[44] G. Prando, P. Carretta, A. Rigamonti, S. Sanna, A. 
Palenzona, M. Putti, and M. Tropeano, Phys. Rev. B 
81, 100508 (2010) 

[45] A. A kbari, P. Thalmeier, and L Eremin, 
|arXiv: 1301. 6497 ^1 (2013) 

[46] L. Pourovskii, V. Vildosola, S. Biermann, A. Georges, 
Europhys. Lett. 84, 37006 (2008) 

[47] A. Kondrat, J. E. Hamann-Borrero, N. Leps, M. Kos- 
mala, O. Schumann, A. Kohler, J. Werner, G. Behr, M. 
Braden, R. Klingeler, B. Biichner, and C. Hess, Eur. 
Phys. J. B 70, 461 (2009) 

[48] J. P. Carlo, Y. J. Uemura, T. Goko, G. J. MacDougall, J. 
A. Rodriguez, W. Yu, G. M. Luke, P. Dai, N. Shannon, 
S. Miyasaka, S. Suzuki, S. Tajima, G. F. Chen, W. Z. 
Hu, J. L. Luo, N. L. Wang, Phys. Rev. Lett. 102, 087001 
(2008) 



